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Abstract 

Threshold photon electroproduction off the proton allows one to measure new electromagnetic 
observables which generalize the usual polarizabilities. There are - a priori - ten "generalized 
polarizabilities", functions of the virtual photon mass. The purpose of this paper is to lay down the 
appropriate formalism to extract these quantities from the photon electroproduction cross sections. 
We also give a first estimate of the generalized polarizabilities in the non relativistic quark model. 



1 Introduction 



In this paper Virtual Compton Scattering (VCS) refers to the reaction 

7* + P^7 + P, (1) 

where 7*, 7 and p are respectively a space-like virtual photon, a real photon and a proton. This 
reaction can be accessed experimentally through photon electroproduction off the proton 

e + p — > e' + 7 + p'. (2) 

In reaction (2) the final photon can be emitted either by the electron or by the proton. The 
first process is described by the Bethe-Heitler (BH) amplitude which is calculable from Quantum 
Electro-Dynamics (QED). The second process is described by the Full Virtual Compton Scattering 
(FVCS) amplitude which, in the one photon exchange approximation, is a linear combination of 
VCS amplitudes. 

Assuming that the one photon exchange approximation is valid, the BH and FVCS amplitudes 
correspond to the graphs shown in Fig.(l). 




(a) (b) (c) 

Figure 1 (a): FVCS amplitude, (b,c): BH amplitudes 



The potential interest of VCS has already been emphasized in the case where the reaction is 
hard enough to allow the use of perturbative Quantum Chromo-Dynamics (QCD) [1]. In particular, 
the interference between the FVCS and BH amplitudes may provide a stringent test of the phases 
predicted by QCD. However, in the hard scattering regime, the experimental difficulty is extreme 
and requires a new type of electron accelerator [2] with high energy (15 -f- 30 GeV), high intensity 
(~ 50 fiA) and high duty cycle (~ 100%). 

In this paper we consider the kinematic regime defined by a center of mass (CM) energy (y/s) 
of the final photon-proton system close to threshold. In particular, we restrict our consideration to 
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values of t/s below pion production threshold. Though the corresponding experiments still require a 
high duty cycle, they are feasible at CEBAF where a preliminary experiment has been scheduled [3]. 
At small Q 2 (minus the virtual photon squared momentum) they can also be done at MAMIB [4]. 

The interest of the region with low energy but arbitrary Q 2 is that, in principle, it allows 
one to measure new electromagnetic observables which generalize the usual magnetic and electric 
polarizabilities. The latter have been the subject of much theoretical interest. They have been studied 
in the context of the constituant quark model [5], the MIT bag model [6], the chiral quark model 
[7], and chiral perturbation theory [8]. For a review on the polarizabilities, see reference [9]. These 
observables are clean probes of the non perturbative structure of the nucleon and are complementary 
to the elastic form factors. By contrast with the latter, the polarizabilities are controlled by the 
excitation spectrum of the QCD Hamiltonian but, because of the threshold condition, the excited 
states contribute virtually. This is an advantage over real resonance production for which meson 
production and rescattering is the dominant process and must be treated explicitly, even though it 
is not the interesting part of the problem in so far as the nucleon structure is concerned. Below 
the pion production threshold, the situation is simpler. In particular it makes sense to neglect the 
strong decay of the excited states, which is the hypothesis of most models of nucleon structure. 
This advantage is somewhat similar to the Deep Inelastic Scattering (DIS) regime for which it 
makes sense to use the parton model because the fragmentation of the partons can be neglected. 

The purpose of this paper is to lay down the appropriate formalism to extract the "generalized 
polarizabilities" (GP) from photon electroproduction cross sections. In short it goes as follows. 

Let q' be the CM energy of the final photon. Then Low's theorem [10] states that, in an 
expansion of the amplitude in powers of q', the terms of order (q') 1 and (q')° are entirely 
determined by the elastic form factors of the proton. As will be shown below there are neither 
ambiguities in the choice of the form factors nor off-shell effects. The terms of order q' have a 
part which is calculable from the form factors and another one which is structure dependent and 
can be parametrized by ten GP's, which are real functions of Q 2 . They are distinguished by the 
spin flip or non spin flip character of the reaction and by the multipolarity of the photons. In order 
to give a term linear in q', the final photon can only be dipole electric or magnetic while, for finite 
Q 2 , the initial photon is only restricted by parity and angular momentum conservation. 
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The paper is organized as follows. In section 2 we specify the notation and kinematics, we 
define the amplitudes and give the expression for the cross section. We shall consider only the 
unpolarized case, the use of polarization being left to future work. In section 2 we also derive the 
general tensor expansion of the VCS amplitude which, in full generality, depends on 12 complex 
scattering coefficients, functions of the photons CM 3-momentum and scattering angle. Section 2 
is not specific to the low energy regime. Section 3 presents the low energy expansion of the VCS 
amplitude. We explain how gauge invariance determines the first two terms of the expansion. In 
section 4 the GP's are defined through the expansion of the non Born part of the VCS amplitude 
over the basis of electric, magnetic and longitudinal vectors. We compute the GP's in a static 
model to illustrate the physical content of these quantities. The scattering coefficients of the non 
Born amplitude are then expressed, in the low energy limit, in terms of the polarizabilities. The 
reason for this intermediate step is that the observables are more naturally expressed in terms of the 
scattering coefficients. In section 5 we tentatively explain how the threshold cross section can be 
analyzed in order to obtain a linear system of equations for the GP's. Section 6 presents an estimate 
of the GP's in the constituent quark model and we conclude in section 7 with a few comments and 
a warning. The various appendices collect the technical material used throughout the paper. 

2 Preliminary definitions 

2.1 Notation and kinematics 

It is possible to develop a covariant formalism [11] but, since we are interested in the threshold 
region, it is natural to work in the CM frame. Therefore, unless explicitly stated, all momenta, 
amplitudes and helicity states will refer to this frame. To simplify the notation the helicity or spin 
labels will be omitted whenever possible. 

The initial electron, final electron and final photon have (4-momentum, helicity): (k, h), (k 1 , h') 
and (q 1 , A') respectively. The initial and final proton have (4-momentum, spin projection): (p, a) 
and (p', a') respectively. In the low energy regime, there is no compelling reason to use the helicity 
to describe the proton spin state and since the (rest frame) spin projection is more usual we stick to 
this choice. We denote as q the 4-momentum of the virtual photon exchanged in the VCS process, 
that is q = k — k' . The spin projections are along q, the direction of q. We note with a roman 
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letter the modulus of a 3-vector, e.g. q = \q\. We choose the coordinate system such that the 
components of the momenta are 

sin 8 sin a cos tp sin ol cos tp 

(3) 

and the (CM) scattering angle of the electron is 8 e = a' — a. The unit vectors of the frame will 
be noted 





"0" 




sin 8 




sin a cos ip 




sin a' cos tp 


q = q 





,q = q 





,k = k 


sin a sin ip 


,k' = k' 


sin a sin tp 




1 




cos 8 




cos a 




cos ck' 



e (3) = q, e (2) = q X g '/ sin 0, e (1) = e (2) X e (3). 



(4) 



In the following we shall neglect the electron mass (m e ) in comparison with its momentum. 
In this approximation, the helicity is conserved at the electromagnetic vertex, which considerably 
simplifies the calculations. Of course one has to be careful in kinematic situations where the only 
surviving mass scale is m e . This occurs when the angle between the electrons, or between the final 
photon and the electron is comparable to m e /k. In practice these kinematic situations are irrelevant 
for the phenomenology of VCS, so we simply make the caveat that the following formulae must not 
be applied in this region. Finally the invariant quantities which appear frequently in the paper are 



Q 2 = -<1 2 , t= (q- q) , s = (p + q) . 

2.2 Cross section 

The S matrix for the (e,e',7) reaction is 

S fi = i(27r)V(p + k-p'-k'-q') T ee \ 



(5) 



(6) 



with the T-matrix a linear combination of the Bethe-Heitler (BH) and full virtual Compton scattering 
(FVCS) amplitudes: 



rpee *y rprSii _|_ rp 



BH , rrFVCS 



(7) 



The unpolarized, invariant cross section is (here the momenta do not necessarily refer to the CM 
frame and m is the mass of the proton) 



da 



(27T)" 



32^/(p.k) 2 -m 2 m 2 e k ° P ° q ° 



d -^4^S\p + k-p'-k'-q')M, 



(8) 
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with 



M = l > IT 

crcr'/i'A' 



-iee 7 



(9) 



The last quantity is invariant under a Lorentz boost and we shall evaluate it in the CM frame. Note 
that eq.(8) supposes that T ee ' 7 is evaluated with states normalized like 

<p\p' >= (2tv)% S{p-p). (10) 

In practice the experiment is performed by detecting the electron in coincidence with the proton, 
and the photon electroproduction event is tagged by a zero missing mass: (k + p — k' — p') 2 = 0. 
The corresponding differential cross section in the laboratory is then 

dai ab = (2tt)~ 5 k' lab s - m 2 ^ 
d\! lah dk\ ah dpf cm 64m k lab s 

2.3 One photon exchange amplitude for photon electroproduction 

In the one photon exchange approximation and in the Lorentz gauge, we have the following 
expressions for the BH and FVCS amplitudes (the relative sign of the BH and FVCS amplitudes 
would change in the case of a positron beam) 

T BH = Z^ £ '- L ^ u{p')V v {p\p)u{p), (12) 

T FVCS = H ^ u(k') lv u(k), (13) 

where e'*(q', A') is the polarization vector of the final photon and e ~ y / 47r/137 is the proton charge. 
There are several choices of form factor decomposition of the vertex Y v which are equivalent for 
on shell particles. We shall prefer Dirac decomposition: 

T»(K',K)=F 1 (u)Y + iF 2 (u) ^ ± u={K'-K)\ ^ 

Fi(0) = 1, F 2 (0) = 1.79, 
because it simplifies the derivation of the low energy expansion of the VCS amplitude. The reason 
is that the Born term of the VCS amplitude satisfies gauge invariance if one uses this decomposition. 
This is not the case for the Gordon decomposition (Gm = F\ + F 2 ) 

Y v {K\K)=G M {u) 1 v -F 2 {u) {K '^ K) . (15) 

Zm 
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From QED perturbation theory one derives the following expression for the (non symmetric) 
leptonic tensor L^ v 



= u(k') 



1 1 

/y ' /y _|_ /y /y ' 

-y.(k' + q') — m e + ie -y.(k — q') — m e + ie 



u 



(k). (16) 



Using the Dirac equation and the anticommutation rules of the Dirac matrices, it is straightforward 
to check that L^ v satisfies gauge invariance, that is 

q'pL'"' = (q v - q' v ) = 0. (17) 

(In the BH process the exchanged photon has momentum q — q'.) 

The hadronic tensor E^ v is actually the object under study. We must therefore avoid 
hypotheses which are too specific to any model. In fact, the only property really needed for 
the phenomenological analysis is that it satisfy gauge invariance, that is 

q'^H^ = H^q v = 0. (18) 

We take this for granted and use it as a constraint on the amplitude. 

To focus our ideas it is nevertheless useful to have a more specific representation of E^ v . For 
this we assume that the interaction of the photon with the hadronic current has the form 

Yy = J df^A^J^r) + ^(r^(f)^(f)}, (19) 

where and are independent of the photon field A^. This form is general enough to encompass 
any reasonable model. The tensor is a contact (or seagull) term which is generally present in 
any model of nucleon structure. It is most often generated by minimal substitution in a quadratic 
kinetic energy term, and thus is necessary for gauge invariance. Its contribution to E^ v is simply 

HTea g ull=<V\S^mV>. (20) 

Using lowest order QED perturbation theory one then gets 



J 'E{ V ') + q' - H n + ie ^ ' 



1 



(21) 



where H n is the strong interaction Hamiltonian. 



2.4 The VCS amplitude 

It is convenient to expand the lepton current in the basis of the polarization vectors 
(e /i (g,A), A = 0, ±1) defined in Appendix B. From the expressions for the helicity spinors 
(Appendix B) one gets, after a little algebra 

u(k', h)Yu(k, h) = J2 n ( h > A ) e "(9> A )> ( 22 ) 

A 

n(h, A) = ( -\e- iX(p a(\h) + S(X, 0)^^27 ) . Q , (23) 
V Q J Vl - e 

+ 2^^/^^7 (k + k') 2 - q 2 

a ( z ) = m ' e = n T^2 o > (24) 

V2 (k + k') + q 2 

where e is the usual polarization parameter in the limit m e = 0. 
We then define the VCS amplitude by 

T FC5 (A',A) = e';(\') H"" e„(\), (25) 

so that we can write 

T FVCS (\') = ^L^fi(M)T FC5 (A',A). (26) 

The VCS amplitude is the main object under study. Its transverse part (A = ±1) evaluated 
along the line q' = q (which implies Q 2 — > due to energy conservation) is the real Compton 
amplitude. (Since the lepton current has been decoupled through eqs.(25, 26) , there is no risk in 
taking the limit Q 2 — > in the VCS amplitude. This is not the case for the FVCS amplitude.) 

2.5 Tensor expansion of the VCS amplitude 

After the gauge degrees of freedom have been eliminated the real photon has two degrees of 
freedom (A' = ±1), while the virtual photon has three (A = 0,±1). Combined with the spin of 
the nucleon this gives 24 amplitudes which are related in pairs by parity invariance. Thus for 
fixed (q, q', 6), the VCS amplitude depends on 12 independent amplitudes. Time reversal gives 
nothing because q / q' - except in the case of real Compton scattering. In the latter case the 
initial photon has 2 degrees of freedom which leads to 8 amplitudes, upon which time reversal 
imposes 2 further relations. 



A convenient way to define the decomposition of T vcs into independent amplitudes is to use 
gauge invariance to write are spatial indexes) 

T vcs = e'{ HijEj, (27) 

with 

£i = £{ - e.qqi/ql- ( 28 ) 
Note that the Lorentz condition has been used to obtain eqs. (27,28). For the real photon there is 
no difference between e' and e* due to the transversality of the (A' = ±1) vectors. 

From eq.(27) we see that T vcs must be scalar with respect to rotations, and bilinear in the 
vectors e '*, e. It is also a 2x2 matrix in proton spin space. Therefore we write 

T VCS = A+ Bia.e(i) (29) 

2=1,3 

with a the Pauli matrices. 

Under a reflection about the scattering plane, a.e (2) does not change while a.e (1), a.e (3) 
both undergo a change of sign. Therefore B 2 is a scalar while B\ and B3 are pseudoscalars. A 
possible choice for the independent scalars and pseudoscalars which are bilinear in (e '*, e) is 

5i = e'*.e, S2 = e'*.q e.q, 53 = e'*.q e.q' (scalars) 

^ ^ / (3°) 

pi = e'*.q e.qx q' , p2 = e'^.qx q' e.q, ps = e'*.qx q' e.q' (pseudoscalars). 

Other combinations either involve e '*.q' , which is zero, or can be written in terms of the or pi. 
Now we consider separately the longitudinal (A = 0) and transverse (A = ±1) amplitudes. 

Longitudinal amplitude: 

Since e^A = 0) = — qQ 2 jq\, we have the relations 

Q 2 ^ 

51 = 2" e .q, 52 = 51, 53 = cos 9 si, 

"° Q2 (3D 

Pi = 0, P2 = T e'*.qxq', p 3 = cos9p 2 . 

% 

From this and the parity requirements we conclude that the most general form of the longitudinal 

amplitude is 

Af- 1 T VCS (X = 0) = a l e'*.q 

r i ( 32 ) 

+ i b[ e'*.qx q' a.e (1) + b\ e'*.qa.e(2) + b\ e'*.qx q' a.e (3) 

where the normalization coefficient Af = ^/^pop^ is introduced for later convenience. 



(33) 



(34) 



Transverse amplitude: 

In this case we have e.q = and e = e, which imply the relations 

->l-k -> n -> '* ^ -* ^' 

S\ = £ .£, S2 = 0, S3 = £ .q E.q , 

pi = e e.g X q\ p 2 = 0, p 3 = e '*.q X e.g'. 
Therefore the expansion of the transverse amplitudes has the form 

A/--1 rpVCS(\ t ->i+ -> | t\ 

+ e'*.qe.qx q' + 6*' e'*.qx q' e.q')a.e(l) 
+ i(b\ e + 6£ e'*.?e.g')ff.e(2) 
+ e'*.qe.qx q' + 63 e'*.qx q' e.q'^a.e (3). 

The 12 scattering coefficients (a 1 , b[ . . .) are complex functions of (q, q', cos 6). Below the 
threshold for pion production the only possible on-energy- shell intermediate states in the RHS of 
eq. (21) are composed of the nucleon and of at least one photon. Their contribution is therefore 
of higher order in olqed an d can be ignored. In this approximation T vcs is hermitian (due to 
the unitary condition) and if we combine this with time reversal invariance we can write, restoring 
the arguments of T vcs ', 

T VCS (q',X',a'; q,X,a)* = T VCS [q\ A', a'; q,\,o) (35) 

where the bar means the time reversed state. From eq.(35) it is then straightforward to show that, 
below pion threshold, all the scattering coefficients are real. 

Finally it may be useful to quote the following relation between the scattering coefficients 
(a 4 , a 1 ') and the usual parameters [12] of low energy (real) Compton scattering (a, /3): 

a 1 ' -> /3/e 2 , a 1 -> - (a + /3cos0)/e 2 when (q = q') -> 0, (36) 

where the limit is taken by first setting q = q' and then letting the common value go to zero. Note 
that (a, (3) are often measured in Gaussian units. In this case one must use e 2 = olqei) ~ 1/137 
in eq.(36) to get our scattering coefficients (which do not depend on the proton charge). 



3 Low energy regime 



The tensor expansion of the previous section is valid at any energy and the scattering coefficients 
are arbitrary functions of (q, q', cos 9). Now we study the low energy regime. More precisely we 
want to define an expansion of E^ v in powers of q'. In this section we review how Low's theorem 
determines the first two terms of this expansion. We propose a rather detailed (re)derivation because 
it helps to define precisely the GP expansion which appears in the next term. 

To define the low energy expansion one needs to separate the non-regular contribution due to 
the nucleon intermediate state in eq.(21). The remainder can then be expanded in partial waves. 
However this separation is not gauge invariant, which makes the analysis of the non nucleonic part 
involved. This is why we add the minimum amount of extra terms to the nucleon contribution so 
as to make it gauge invariant. Those extra terms are of course subtracted from the non nucleonic 
contribution so that the total E^ v is unchanged. 

3.1 Preliminary lemma 

We shall say that an amplitude is regular with respect to a generic variable q if, in the vicinity 
of q^ = 0, it can be expanded in powers of the 4-momentum q^ (rather than in powers of q and 
q only). Thus, if T is regular we can write 

T = T(°) + T,? V + T$qi>q v + ■■■ (37) 

where the coefficients (T^°\ T^ 2 \ . . .) are independent of q^ . For instance, a.q is regular 

while a.q/q is not, even though it is finite at q = 0. Note that what matters is that the amplitude 
can be written in the form (37). The fact that the energy component qo can be eliminated in favor 
of q is irrelevant at this stage. 

We now use the notion of a regular amplitude to establish the following lemma. 

Lemma 

Let us assume that an amplitude, Tp, for the production of a real photon is regular with respect 
to its momentum q'^ and satisfies a gauge invariance condition : q'^T* 1 = 0. Then is of order 
q' at least. 
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Proof 

For a real photon we have q' = q'. Therefore we can write 

T» = a> 1 + 6^^ + 0(q' 2 ). (38) 
where (q' 2 ) means terms of order q' 2 or more. From the gauge invariance condition we have 

= q'y + O (q' 2 ) = q'a + q\a { + O (q' 2 ) . (39) 
Since the coefficients are constant we can average the above equation over q ', which leads to 

q' a ° + 0(q' 2 ) = 0. (40) 
We can also multiply eq. (39) by q' and average again. This gives 

qV + 0(q' 2 ) =0, i = 1,2,3. (41) 
Since is independent of q', eqs.(40, 41) imply that = 0, which proves the lemma. 

3.2 Born and non Born amplitudes 

To apply the lemma we must first separate the non regular part of E^ v . This can only come 
from the nucleon intermediate state in eq.(21). The higher energy states start at *Js = m + m T , 
which guaranties that the energy denominators do not vanish at the soft photon point. So at least 
term by term the sum over all the possible states but the nucleon is regular. The seagull term, 
if present, is also regular since it is the matrix element of a local operator. Therefore a natural 
starting point is the decomposition 

= Hg + BT + HZ aguW (42) 

with HpJ' the nucleon contribution, and R' 1 " the contribution of the higher energy states. From 
eq.(21) we have 

*jr =u( P ')r»( P ',p>) J[ P 'J r(p', P ).( P ) 

+ ^' f » 2P,(^-f,) r(f -' ),M ' 
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where we use the notation Il(K ) = u(K )u(K ) and 

r^= £ <p'\J' 1 (o)\x,p , > 1 ^^<x > pV(o)Ip> 

jfev 2E x (p')(p + q'-E x (p')) 

v yv v yy (44) 

2E X (P J(po -q'-E x (p)) 

with P' = p' + q', P = p — q' , Exij^j = \J K 2 + m^-, mj being the mass of the excited 
state. Note that the sum in eq.44 includes the nucleon-antinucleon excitation, or pair term. In 
eqs. (43,44) the intermediate states are on mass shell. So at this point the choice for the vertex T* 1 
is still arbitrary. We denote as P' 11 = (VP' 2 + m 2 , P'^ and P? = {^J P 2 + m 2 , P^ the on mass 
shell 4-momenta of the intermediate nucleon in the direct and crossed terms respectively. 
Explicit calculation shows that 

q'pHg^O, #£W0. (45) 
The reason is twofold. First, the intermediate 4-momenta are off the energy shell, that is 

P'fi + Pfi + l'fi, Pfi^Pfi- q'p- (46) 

Second, the nucleon state alone is not enough to have closure in the space of Dirac indices. This 
would prevent gauge invariance even for a point-like particle. The antinucleon contribution is 
necessary for completeness. To restore gauge invariance we proceed in two steps. 
First, we write 

H N = H N + (H N - H N ) (47) 

where is obtained from by replacing, in the vertex T only, the 4-vectors P^ and P^ by 
(p' + q')^ and (p — q')^ respectively. We have 

H% - H% = u{p') 



+ 



u(p) 



2P' Q {p' Q + q' - P' Q ) 2P (PO -q'-Po) 

d"" = t»(p',p') n(p') r (p', P ) - T»(p',p' + q ') n(p') + q \ T ) m 
c^ = t v {jp\p)u{p) r"(p , v )-r{p\p-q') n(p) r"(p-g' >P ) 

Since 

Pi- (p' + g')^ = %.°)(^o-Po-q'). P »- (p-q^^S^^iPo-Po + q), (49) 
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the energy denominators in (jIn — Hn^J are exactly balanced at q' = by the differences in D^ v 
and C^ v . Therefore (jIn — Hn^J is not only finite but regular with respect to q'. 
Second, we have u(K)u(K) = 7. if + m, and from the identities 



7.P' + m _ 7 .(p' + q') + m ~j P' Q + 7.P 



(50) 



(51) 



2P'M + q' - i*) (pi + g') 2 _ m 2 2P ( J(p( ) +q' + i*) 

-y.P + m _-y.(p — q') + m -jqPq + 7.P — m 
2P (po - q' - Po) ~ (p_ g ') 2 _ TO 2 " 2P (p -q' + Po) ' 

we see that we can write 

H N = H B -H„, (52) 

with the Born term defined by 



h% = «(j/)r"(pV + g ') jf' + ' 2 ' ) + T ! r"( P ' + g ', P )u( P ) 



(p — q) — m 



(53) 



while the antinucleon contribution , which comes from the second terms in eqs.(50,51), is 
manifestly regular. In summary, we have 

H N = H B -H W + (h n - H N ) =H B + (H n - H B ) (54) 

and what we have just shown is that (iijv — Hb) is regular with respect to q'. Obviously this is 
also true with respect to the momentum of the virtual photon. 

3.3 Gauge invariance of the Born term 

The last point is to establish the gauge invariance of Hb as defined by eq.(53). Until now we 
have not used the explicit form of the vertex T. The only restriction is that the same form must be 
used in the Born term and in the non Born term, because it is essential to establish the regularity 
of the latter. Now we use the vertex decomposition of eq.(14) from which we find 



qfl H B =u(p) -y.q ^ - T (p+q)+T (p,p-q) -y.q 



u(p). (55) 
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Using the anticommutation relation of the Dirac matrices and the Dirac equation this can be 
transformed to 

= u{p') [T v {p' + q',p) - T v {p',p- q')]u{p) = 0. (56) 

Note that the cancellation between the two terms is due to the fact that the vertex defined in eq.(14) 
depends only on the difference of its arguments. 

A similar derivation leads to H B ' q v = 0. Apart from the limit Q 2 = 0, this is not useful as far 
as the low energy theorem is concerned, but it is necessary for the forthcoming developments. 

3.4 Low's theorem 

We now have the elements to derive Low's theorem. We define a new decomposition of the 
amplitude 

H = H B + H NB , H NB = R + (H N - H B ) + H seagu u, (57) 

where Hnb is regular - since both R, H seagu u and (Hn — Hb) are regular. Moreover, the gauge 
invariance of H and Hb implies that H^b is itself gauge invariant. Therefore the lemma tells 
us that Hnb is of order q'. Since, for finite Q 2 , Hb is of order 1/q', the first two terms of the 
expansion of H in power of q' are determined by the known Born term. 

4 Generalized polarizabilities 

In this section we deal with the next term, that is the term of order q', in the expansion of 
the amplitude. Part of this term comes from the Born amplitude, but this part is trivial and can 
be computed exactly. So we set it aside for the moment and consider, in this section, only the 
non-Born part of the amplitude. 

4.1 Partial wave expansion of the non-Born VCS amplitude. 

To expand Hjf B we use the complete basis of 4-vectors {V^^pLM, q), p = 0, . . . , 3) given in 
Appendix C. At this stage we need not assume that the final photon is real, so we keep all possible 
values of p' for symmetry reasons. We thus write the expansion in the form 

K V B {q'<r',q<r) = ^ X V (p' 'L' M',q') H^' '' pLM '(q» V, q<r) g pp V»*(pLM,q), (58) 



14 



where the sum runs over all possible values of (p'L'M', pLM) and where we have restored the 
full set of arguments of Hnb- Using the orthogonality of the basis vectors we can invert eq.(58) 
and get the expression for the multipoles 

< tf '^(qV',q^) = (47TJV)- 1 J dqdq'V^p'L'M', q') H^ B (q'*', qa) V u (pLM, q). (59) 

From the relations (see Appendix C) 

q ll V' l (OLM,q) = q Y&(q), q^SLM, q) = -qY^q), q^V^^pLM, q) = 0, p = 1,2 (60) 

we see that gauge invariance imposes the following constraints 

0, 

(61) 



,rr3L'M',pLM , ttOL' M' ,pLM _ n 



n p'L'M',3LM n p'L'M',0LM _ n 

These constraints are satisfied by eliminating the multipoles with p = 3 in favor of the ones with 
p = 0. If we define the new basis vectors 

W^{pLM, q) = V^ipLM, q), p = 1, 2 

on (62) 
W^iOLM, q) = V^^OLM, q) + —V^^LM, q), 

q 

which satisfy (q fl W' 1 (pLM, q) = 0, p = 0,1,2), we can write the manifestly gauge invariant 
expansion 

H^ B (q'a',qa) = g, p W»(p [ V M',q') H&* M ' ' > pLM (qV, qa) g pp W v \pLM,q), (63) 

where the sum is now restricted to p, p' = 0,1,2 (and, of course, still runs over all possible values 
of L, V , M, M'). The usual nomenclature for these multipoles is : p = 0: longitudinal or charge 
type, p = 1: magnetic type, p = 2: electric type. 

From their definition, eq.(59), and the tensor character of ' , it is easy to show that the 
multipoles have the following property 

H P N L B M '' pLM (q'a',qa) = ^ < L'M'\1l\L'l^ X ^'1^1^ > 

a,a',M,M' " (64) 

< LM\TZ\LM X - a\R\ -a > H P N L B W ' ' pLW (q'a' , q a) , 



15 



where 1Z is the rotation operator. By integrating both sides of eq.(64) over the rotation group and 
by using the well known formula for the integrals of rotations matrices [13] one can factorize the 
dependence on (M'a'Ma) into Clebsch Gordan coefficients. If we define the reduced multipoles by 

rr(p'L',pL)S ( , \ _ 1 \- 

a'aM'M (65) 
(-) 2 < - - a , -a\Ss >< L M ,L - M\Ss > H NB ,H (qtr,qtr), 

then 

rrp'L'M',pLM til nrr \ _ i s\W+L+M 

h nb (q a >q°v - {-) 

Ss (66) 
<\- *\ l -*\Ss >< L'M',L- M\Ss > H^^ pL)S {q\q). 

This completes the angular analysis. By combining eqs.(63, 66) one can write Hnb in terms of the 

reduced multipoles defined by eqs.(59, 65) - all the angular and spin dependence being collected 

in known functions. From eqs.(59, 65) one has the following selection rules 

S = 0<xl, \L' -S\<L<L' + S, (-) P ' +L ' = (-) P+L (p',p = 0,1,2), (67) 

and using the time reversal properties of the hadronic tensor, it is straightforward to show that, 
below pion production threshold, the reduced multipoles are real. 

4.2 Low energy behavior and Generalized Polarizabilities 

As we want to define threshold quantities at finite q, we must consider the limit q' — > at fixed 
arbitrary q, which we call the low energy VCS limit. At q' = 0, the virtual photon has energy 
qo = m — *sjm 2 + q 2 . Therefore if we let q go to zero, the initial photon becomes soft but it 
remains space-like because \qo\ ~ q 2 /2m < q. As a consequence the limit q — > of VCS may not 
coincide with the limit along the real photon line, that is (q = q') — > which is appropriate to real 
Compton scattering. For instance if we consider the Born term to leading order in q', we get 

T B ~ -^e'*.e (VCS limit) (68) 

q 

Tb ~ e'*.e (Compton limit) (69) 

These two equations are formally identical since in the Compton limit go/q' = 1, but 
experimentally this makes a difference since along the VCS limit \qo\ ~ q 2 /2m. As far as the 
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Born term is concerned, this does not matter since we know how to compute it exactly. The 
problem becomes serious when one tries to define the polarizabilities. Since these are experimental 
quantities their limit when q — > will actually correspond to qo ~ q 2 /2m and therefore may differ 
from the same quantities measured in real Compton scattering. To resolve this problem we need 
to analyse the low energy behavior of the non Born multipoles. 

To this end we use the fact that H^ B has an expansion in powers of the 4- vectors q and q'. 
For arbitrary (qo, q' Q ) elementary harmonic analysis shows that 

X =4r~l I dq'dqY l \q")H NB {q', q)Y l (q) = G(qi qo) + 0(q)0(q') (70) 

q " q J 
where G has the form 

G(qi qo) = a(l', I) + b(l' , l)q + b\l\ l)q' Q + ... (71) 

When we substitute the actual value of qo and q' we find that the limit as (q, q' — > 0) is independent 
of the path since we have in any case (qo,q'o — > 0) so that we always end with G = a(l',l). Of 
course, this assumes that the expansion of G converges fast enough but this can be seen as a 
consequence of the fact that Hnb is regular. 

From the above analysis we see that quantities of the form X have the desired property, that 
is they have the same limit as (q, q' — > 0) along the VCS path and along the real Compton path. If 
we go back to the non Born multipoles ,pL ^ S , there will be no problem when (p' , p = or 1) 

because the rank of the spherical harmonics of q or q' are just equal to the total angular momentum 
of the photons L or V . So we can define the GP at arbitrary q as 

p',p = 0orl. (72) 

q'=0 

If the final photon is of electric type (p' = 2) we can use the fact that we are only interested 
in the leading behavior of the multipoles in q'. This allow us to relate them, for any p, 
to the GPs already defined. Using the definitions of Appendix C and taking the appropriate 
linear combination of the electric and longitudinal vectors one derives the following relation ( 

v^l'm'a') = (o,e L , M >(q')) ) 

W(<7') = \j J ^YV l {q<)q< + yj 2 JL^ly^^). (73) 



p( P 'L', P L)S^ 



q' L 'q L NB 



i^/^V, q) 
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When averaging over q' to compute the multipole, the second term in the previous equation gives 
a contribution of order q' L +1 . The first term can be transformed using gauge invariance, that is 



and one gets 



^liTjiu % ttOv ttOv (HA \ 

1 H NB = -J H NB = H NB ( 74 ) 



H^" L)S W, q) = -^^lf W(q\ q) + 0(q' i '+ 1 ) (75) 



which is known as the Siegert [14] relation in other contexts. If (p = or 1) we get back to the 
cases already studied. So we are left with (p = 2,p' = 0,1,2) but we can still use the Siegert 
relation in the case p' = 2. So we only need to study the case (p' = or 1, p = 2). Using eq.(73) 
in the case of the virtual photon, and gauge invariance of Hnb 

Kb? = q jK°B (76) 

we get the relation 

(p'L',2L)S _ ITTTqo (p'L',0L)S j lL + l ~(p>L',L)S 

NB L ~q NB ~V L NB 

where Hnb is defined from 

H(r'L'M',LM) = J d ^ Vfi{p , VM ^ *, )H g B $) ] 1 (V8) 

Note that Hnb is neither electric nor longitudinal. 

From eq.(77) we see that, when p = 2, we cannot define a GP which has the desired photon limit 
when q — > 0. Indeed from the results for the case (//, p = or 1) we know that when (q, q' — > 0) 

^' ,0i)5 (q', q) ~ q'^q 1 + 0{^ L '^)0{^) (79) 

and obviously from eq.(78) 

H { ^'' L)S ~ q'V +1 + 0(q' i,+1 )0(q i+2 ) (80) 

with a path independent limit. In the real Compton case we have qo = q, so that the first term in 
eq.(77) is leading, while in the VCS case \qo\ ~ q 2 /2m and the two terms are of the same order in 
q. To preserve the real photon limit we therefore define the following mixed GPs 



p( P 'L',L)S 



q iL' q L+l^NB 



tffcf' i)s (<.',q) 



(81) 

q'=0 



18 



Collecting the results of the above discussion, we can now summarize the low q' behavior of all 
the multipoles at arbitrary q: 



H% L B > PL)S W, q) = q' V W(q) + 0(q' i '+ 1 ) (p = 0, 1) 



(82) 



^L+±kpilL', L)S {q) + q ^TTp(lX' ! X)5 (q) 



\V L 

q q 



+ 0(q' i,+1 ) 



'^V, q) = -^'q V ^P^'^M + 0(q' i '+ 1 ) (p = 0, 1) 

^'' 2i)5 (q',q)_= 

^ + T| p (0X' ! i )5 (q) + q y^±lp(0X',X)5 (q) 



,L> L L' + l 

q q 



+ 0(q' i,+1 ) 



(83) 



(84) 



(85) 



We have omitted the case p' = since it is not needed for a real final photon. If we now specialize 
to the multipoles which are linear in q', that is V = 1, we can use the selection rules, eqs.(67), to 
enumerate the possible quantum numbers of the multipoles. They are shown in Table 1 and we see 
that, to describe the low energy behavior of the VCS amplitude, we need 10 GPs which are 



p(ll,00)l p(ll,02)l p(ll,H)0 p(n,n)i p(H,2)l 



(86) 



p(01,01)0 p(01,01)l p(01,12)l p(01,l)0 p(01,l)l 



(87) 



Final photon 


Initial photon 


Spin flip type 


P' = l 


p = 


L = 0,2 


S = 1 


p = 2 


L = 2 


S = 1 


p=l 


L = 1 


S = 0,1 


p' = 2 


p = 


L = 1 


S = 0,1 


p = 2 


L = 1 


S = 0,1 


p=l 


L = 2 


5 = 1 



Table 1 Allowed quantum numbers of the non Born multipoles, to order q'. 
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4.3 Expressions in a static model 

To get some insight into the physical content of the GPs defined above, it is useful to have their 
expression in some simplified situation. For this we consider the limit of a very heavy nucleon 
so that all velocity effects, which go like q/m or q' /m, can be neglected. For simplicity we also 
ignore possible seagull terms. 

To compute TV 1 " from eq.(44) we need matrix elements of the form 

< X'K'IJ^O^X K >= e i ( # '- # )- f < X'K'IJ^IX K > . (88) 

The state \X K > is obtained by boosting the state |X0 > to velocity K/Kq. Thus the difference 
between \X K > and |X0 > generates only a velocity effect (recoil) that we neglect. 

The translationally invariant state |X0 > is proportional to the intrinsic state \X > normalized 
to unity: 

\X0>=x\X>, <X\X>=1. (89) 
From the normalization convention of eq.(10) we get, in the zero velocity limit 

x 2 = (27r) 3 2mx s(o^j = 2m x V, (90) 
where V is the normalization volume. Therefore we can write 



<X'K'\Ji i (0)\XK > = ^ J dfe i ( g '-^)- f <X'K'\J' 1 (r)\XK> 



(91) 



which leads to (c.f. eq.44) 



<N\ j df' e-^'- f >(f')|X >< X\ J dre i ^r(r)\N> 



2m J J m- mx + q' 



+ <N\ f dfe ilF J"(r)\X >< X\ I dr l 'e~^ ' y ' >(r')|W > 

J ' " J m-m_ 

where recoil effects in the energy denominators have been neglected. 



m — mx — q' 



Here we consider only two significant examples, which show the connection with the classical 
polarizabilities. Starting from the definition, eqs.(59, 65, 72), and using eq.(92) we find : 

P (oi,oi)o (q) = J|^[< N\d(0)\X >< X\d( q )\N > 
x 

+ < N\d(q)\X >< X\d(0)\N >} _ l — + £p( 01 - 01 ) , (93) 

d(q)= [dr 3 -^Ar) Z , 
J V 



m — mx 
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and 

p(11 ' n)0(q) = 7f5> N ^ x ><: X ^ N > 

+ < N\fi(q)\X >< X\fi(0)\N >] 1 + 8P( n ' n )°, (94) 

m — mj 

2 J qr L J z 

Here is a spherical Bessel function and SP denotes the part of the polarizability which comes 
from (ifjv _ Hb) in the definition of Hj^b, eq.(57). In this static model its effect is simply to 
cancel the anti-nucleon pair excitation from the sum over states in R. If we agree to not include 
the anti-nucleon excitation in R, then we can set SP = 0. 

When q — > 0, [3ji(qr)/qr] — > 1 so that, up to a constant factor, we recognize in eqs.(93, 94) 
the usual electric (a) and magnetic (J3) polarizabilities of the nucleon [12]. Explicitly: 

P (01 ' 01) °(0) = -yf <*/e\ P^ n )°(0) = -yf /3/e 2 . (95) 

In the examples we have shown, it is clear that the virtuality of the photon means that we 
actually measure a transition form factor. In so far as a given GP is dominated by a single excited 
state, measuring this quantity amounts to measuring the corresponding transition form factor without 
really producing the state. This is likely to be the case for the A(1232) which dominates p( 11 ' 11 ) 
because of the absence of radial dependence in the soft photon vertex. 

4.4 Relation between scattering coefficients and 
generalized polarizabilities 

The observables are most easily calculated in terms of the scattering coefficients defined in 
section 2.5. In this section we relate these coefficients to the GPs by a method analogous to the 
one used in the analysis of pion electroproduction reactions [15]. We define 

1 B — £ fi n B £j/ > 1 NB — £ ti n NB £l, -> vyt>J 

and we expand Tj^s m trie same way as the full amplitude, eqs(32, 34), with coefficients noted 
(a l NB >•••)• If necessary, the corresponding coefficients for the Born amplitude can be obtained 
from eq.(53). Here we consider only the non Born part. 
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Using expansions (63, 66) we write 

p',p=0,l,2 L',M',LM,S 
g plpl E , ;{\')W»{p'L'M)g pp e fl {\)W»(pLM,q) (97) 



( _ ) L +CT ' +i+ M < I _ ^ 1^1^ ><; L » Mfj ^ _ M \ Ss > 

2 2 



With our choice of polarization vectors (see Appendix B), we get (e s = Q 2 /q q) 

e^(0) W*(pLM, q) = 



e s Y M *(q) ,p = 
,p = l,2 

-e(±l).Mj$(q) ,p = l (98) 
£ „(±1) W*(pLM, q) = \ _£(±i).£ , p = 2 

,p = 

and similar relations for the real photon. Now we have 

M M (q) = C L L q Y^(q), £ M (q) = -iC L qX L q Y^(q), (99) 

where L q = —iqxV q and Cl = 1/ \/L{L + 1). If we substitute eqs.(98, 99) into eq.(97) and take 
into account angular momentum and parity selection rules, we get for the longitudinal amplitude: 

£ c v e».l t (h^ 1 '- 1)1 mL^ 1 + j/ft'^ m!^ 1 ' 1 ) (iQo) 

i',5 

and for the transverse ones 
(47TJVT 1 T^(A',A = ±1) = 

- x L ql r.L^^^^'^'-^Mg?'- 1 ) 1 + C £ , +1 ^'' 1£ ' +1 ) 1 Mif^' +1 ) 1 ) 



L'S 



i',i 



(101) 

where Mi ( ^ ; is an operator in spin space defined by 



i ,i 



<a'\M^ S \a> = (-y + °' +L+M < 1 --a\ 1 -a\Ss><L'M\L-M\Ss>Y£,(q')Yit(q). 
MM' 

(102) 
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After a little algebra the last quantity, considered as a matrix in the spin space of the nucleon, 
can be written in the form 



(103) 



where P £ is a Legendre polynomial. After substitution of eq.(103) in eqs.(100, 101) it then remains 
to evaluate the action of L q and L q i and to project on the tensors used in expansions (32, 34) to 
get the relation between the non Born scattering coefficients (a l NB . . .) and the reduced multipoles 

^ Vi)5 (q',q). 

The above method can be used to get the general angular decomposition of the non Born 
scattering coefficients in terms of the GPs. The resulting expressions are quite long and not really 
useful in our context since we need only the first term in an expansion in powers of q'. This 
amounts to retaining only the terms with L' = 1 which leads to the following relations: 

Longitudinal coefficients : 



J 

a NB 


= A' + 0(q' 2 ), 


a l = 


NB,1 


= /3<q' + 0(q' 2 ), 


P[ = 


NB,2 


= /3k' + 0(q' 2 ), 


A = 


NB,3 


= /3k' + 0(q' 2 ), 


A = 



^qP (01 ' 01) °, 

£sVS ^os^W+q^ 11 ' 02 )^) - ^qP* 01 ' 01 ) 1 
p(ll,00)l + q 2p(ll,02)l /v ^ _ ^ cog g qp (01,01) 

p(ll,00)l _ ^ q 2 p (ll,02)iy 



2sin0 
2 sin 6 



(104) 
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Transverse coefficients 



t t I i r\( '2 

a NB = a q + U(q 



a NB = a q + 0(q 
^ B)1 = /3k' + 0(q 
^i^i'q' + ^q' 

fe jvs,2 = /^k' + o(q' 
fe jvs,2 = /^V + °(q 



12 



a 



a 



^ (q cos ^P( 11 ' 11 )° + 25oP (01 ' 01) ° + V6 q 2 P^ 01 ' 1 ^ ) , 
/3 D (11,11)0 

: -V8 qP ' 

2 ~p(01,01)l + ^ q 2 ^3 p(01,l)l + p(01,12)l N 
4 sin V \ ' / 



4sin0 
3 



, fi\ = sin 2 (9/3f, 



4sin 2 



q p(".")l _ ^ 5oq P (11 '° 2)1 + 2gb cos flpt 01 ' 01 ) 1 



+ V^q'cosflf-Pt 01 ' 12 ) 1 + V^P* 01 ' 1 ) 1 ) - VsT^q 3 ^ 11 ' 2 ) 1 



4 sin 2 



cos 9 qP( 11 ' 11 ) 1 - V372 5bq cos B P^ 1 + 25bi ,(01 '° 1)1 



+ V2 q (-pC 1 - 12 ) 1 + ^P^ 01 ' 1 ) 1 ) - y/5/2 cos q 3 ^ 11 ' 2 ) 1 

(105) 



Note that, to shorten the expressions, the dependence of the GPs on q has been omitted in 
the above equations. 

5 Analysis of threshold experiments 

The aim of this section is to investigate how the (e,e',7) cross section can be analyzed in 
order to extract the GPs. For simplicity we consider non polarized experiments. We shall see 
that, in this case, only 4 independent combinations of GPs can be measured, but given that this 
field is just beginning, the analysis presented hereafter is sufficient to interpret the first generation 
of experiments. The case where polarization is measured is a straightforward extension of the 
material presented in this section. 

5.1 Parametrization of threshold cross sections 



We recall that the unpolarized cross section is 

dai ab (2tt)" 5 k' lab s - m 2 



dk lab dk lab d P 



cm 



64m kiab 



M, 



(106) 
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where M is invariant under a Lorentz boost. We choose as independent variables (q', q, e, 8, <p) 
and denote with a tilde the value of the dependent variables at q' = ( e.g. qo = m— \/m 2 + q 2 ). 
We assume that at fixed (q, e, 0, <p) the experiment determines in the form 

yu ex P yu ex P 

where the coefficients are functions of (q, e, 5, <p). 
We write 

r ee' 7 = T Bff + T FVCS + (108) 



and define 



4 



1 2 A/f-^ 7 Ayfi- 57 

l^ rB H + T |VC5 = ^_ + ^_ + ^ + 0(q »). (109) 



q q 



The superscript LET reminds us that it is the part of the cross section which depends only on the 
low energy theorem, that is the coefficients (A4f ET , i = —2,-1,0) are calculable functions of 
(q, e, 8, <p) once the form factors of the proton are known. The corresponding expressions are very 
long and, since it is relatively trivial to compute T BH + T^ vcs at any kinematic point, it is easier 
to evaluate (Mf ET , i = —2,-1,0) by a numerical procedure. 

Obviously the coefficient Mo exp is the sum of Mq^ 7 and of the interference between T^ cs 

B 

we define 

U = [q' (T BH + T EVCS ) 
V = q'" 1 T NB 



(which is of order q') and the singular part of T BH + Tg VCS (which is of order 1/q'). Therefore 



q ' =0 (110) 

(-1 rrFVCS 



q'=0 

The output of the experiment is then, first, a test of the low energy theorem in the form 

M L _ ET = A4 e _ x 2 p , M L _ ET = M e *f, (111) 
and, second, a linear system for the GPs 

I ^ u\y + UV \ = _ M LET m (112) 

spins 
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To work out the LHS of the previous equation, we need the singular part of T BH + Tg VCS . 
Starting from the definition of the amplitudes (see sections 2, 3), it is a simple exercise to derive 
the following expression 

2m 



U 



e 3 q 



2m 



{ e^27^(g 2 ) 



^/T^7 QY m + £(q) I Q 
Gm(Q 2 ^J [(zhy/l — ecos <p — i\J\ + e sin ip^J a.e (1) 
+ + e cos (y5 + 2^\/l - e sin <p^ a.e (2)] 



(113) 



where 



K 



q k 1 

+ 



k 



(114) 



q'=0 



\ p.g' k' .q' k.q' 
E(q) = Vm 2 + q 2 , and G £ (Q 2 ) = F^Q 2 ) - F 2 (Q 2 )Q 2 /4m 2 . 

Next we need the coefficient of q' in T^ cs . For this one has just to combine eqs.(26, 32, 
34, 104, 105) to get 



V=^-^E(q) Y, n(h,\)A(\), 



Q 



A=0,±1 



with 



and 



A(0) 



a l e'\q 



+ i\(3 l 1 e'*.qx q' a.e (1) + /3 l 2 e'*.qa.e(2) + pl e'*.qx q' a.e (3) 
A(±l) = a e'*.e + a" e'*.qe.q 



(115) 



(116) 



(117) 



+ i(Pl e'*.qe.qx q' + fil'e'+.qx q ' e.q ')a.e (1) 

+ i(& e'*.e + $ e'\qe.q')a.e (2) 

+ i{f3\e'*.qe.qx q' + 0$ e 1 * .q X q' e.q') a.e (3). 
(The dependence of the polarization vectors on the photons' helicities has been omitted to shorten 

the last two formulae.) 

The rest of the calculation is straightforward and leads to 

4me 6 q 



AC P - M L Q ET 



2£(q) 

" Q 2 (l -e)\l m+E(q) 

{sin d(w" sin 6 — u/kr cos (p cos 6) [ePii(q) — -Prr(q)] 

— (a;" sin 9 cos ip — u/kj' cos 6) -\/2e(l + e)PiT(q) 

- (a;" sin 6 cos 6 cos <p - u/k T (l - cos 2 <p sin 2 0))y/2e(l + e)P' LT (q)}, 



(118) 
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where the various kinematic coefficients are 



u> 



p 1 .q 1 k.q 1 
k T = Q 



UJ 



q'=0 



k 1 .q' k.q' 



q'=0 



2(1-0 ' 



u>" = ujq — w'vk' 



(119) 



and finally 

3 
2 



^TT(q) = \Gm(Q 2 ) [2$oP (01 ' 01)1 (q) + q 2 (\/3 i* 01 ' 1 ) 1 + P^ 1 ) 



3 mq y 

2 g"' 



(120) 



As announced, only 4 independent combinations of GPs, that is Pll, Ptt, Plt, P'lt> can 
be measured in an unpolarized experiment. Moreover, to separate Pll and Ptt a measurement 
at 2 values of e is necessary. Since, up to a known factor, Pll is exactly the generalization 
of the usual electric polarizability (see eq.93), an experimental effort in this direction seems 
valuable. As is clear from eq.(l 18), the angular distribution depends only on 3 unknown coefficients: 
ePll + Ptt, Plt, P'lt- The interference between the BH and VCS processes is entirely controlled 
by the kinematic coefficients a>, u/, uj" which have a known ( though non-linear) dependence on 
the angles and energies. Therefore, at least in principle, only a few (a priori 3) angular set-ups are 
necessary to extract ePll + Ptt, Plt, P'lt- 

6 Generalized polarizabilities in a quark model 

The aim of this section is to provide a first estimate of the GPs in a reasonably successful 
model. It is defined by the non relativistic quark model (NRQM) Hamiltonian 

2 



Hnrqm = ^ 



2m 



+ V(r), 



(121) 



where m q is the constituent quark mass and (f n , —iS7 n ) are respectively the position and momentum 
of the quark n. We introduce the coupling to the photon by imposing the condition that the 



27 



Schroedinger equation 



d 



(122) 



should be invariant under the gauge transformation 



r2, r3,i) — > exp 



-iy~]qnv(r n ,t) 



*(ri,r2,r 3 ,t), A^t) -> A^t) + d^v, (123) 



where q n is the charge of quark n in units of e and v an arbitrary function of (f, t). This is achieved 
by the following minimal substitution 
d d 



dt dt 



+ * E qnAo(r n , t), - iq n A(f n ,t). 



(124) 



In order to obtain the correct spin dependence of the interaction we write 

&n-Vn ) ( &n-Vn 



hnrqm = y] 



2m 

71 * 



+ v(f B ) = 



2m 



+ V(f m ), (125) 



before the substitution defined by eq.(124). This trick produces the same result as the non relativistic 
reduction of the Dirac current. A little algebra then leads to the following interaction between the 
photon and the quarks (at time t = 0) 



v ii = Yl ^ qn 

n 



A (f n ) + - (D(r n ) + i* n .B(r n ) + 2A(r n ).Vr 

Zm q 



Zm q 



where 



D = \7.A, B = \7xA. 



This can be written in the form 



V 7 = J drA^J^r) + lA^A^S^r) 



(126) 



(127) 



(128) 



where integration by part has been used to eliminate D and B in favor of A and 



2m„ 



S(f- r n ) I - a n x V m 



+ &n x Vn S(f- r n ) 



(129) 



<^(f) = £ -f^j), = (1,2,3), 



s o» = sll0 = Q _ 
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To simplify this first estimate of the GPs we again neglect, as in section 4.3, contributions which 
vanish with the velocity of the initial or final nucleon. The calculation of these recoil corrections 
is postponed to future work. To perform the calculation of the GPs we need the multipoles of 

Hnb = R + H seagu u + H N -H B (130) 

with the quantum numbers displayed in table (1). According to eq.54 this is also 

Hnb = (R- H w ) + H seagull + H N - H N . (131) 

As explained in Sec. 3. 2, R includes all possible states but the nucleon. In particular it contains the 
nucleon- antinucleon excitation which is exactly Hj^. Therefore the term (i? — Hj^) contains only 
the sum over the resonances and this is what we compute with the quark model. This simplification 
is, of course, due to the fact that we have moved Hjj into the Born term so as to make it gauge 
invariant. Moreover in the non relativistic quark model, the vertex T(K ', K) involves form factors 
which depend only on (jt — K'^j . Therefore the form factors play no role in the difference 
Hn — Hn and when we average the latter which a spherical harmonic of q 1 with rank 1, as is 
required for the evaluation of the polarizabilities, dimensional analysis tells us that the result is of 
order q'/m. This is a velocity effect that we neglect. So within our approximation, only the sum 
over resonances (R — Hj^) and the seagull term contribute. The relevant formulae are collected 
in Appendix D. 

Though we have neglected all recoil terms to establish these equations, we have kept the q 
dependence of the energy denominators of the crossed term. This allows a cheap estimate of the 
size of the recoil corrections as q increases, by comparing with the results obtained when q is set 
to zero. We find that up to q = 0.5 GeV, the effect is always smaller than 5%, and can reach 
20% at q = 1 GeV. Since this represents a second order correction, this is a warning not to trust 
our estimates beyond a few hundred MeV. In the figures we have plotted the results obtained with 
q = in the energy denominators. 

Our results are shown in fig. (2). In this model only seven GPs are non zero, which is due 
to the lack of deformation and to the non relativistic approximations. When the final photon is 
magnetic, the corresponding operator is just the magnetic moment, which has no radial dependence. 
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As a consequence, only the A(1232) contributes to the sum over intermediate states and the non 
relativistic charge density, which has no spin dependence, cannot connect the nucleon and the 
A(1232). This implies pt 11 - 00 ) 1 = pC 11 - 02 ) 1 = o. 

Next, p( 11 ' 2 ) 1 involves a matrix element of the form 

J dqy 2 ^(q) J dre*" < A\J(f)\N > . (132) 

0.00 i . , . , . , . , . 1 0.10 | . , . , . , . , . 1 
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Figure 2 Generalized polarizabilities in the NRQM. 
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Due to angular momentum conservation, this is zero in the absence of deformation in the nucleon 
or delta states, which is an hypothesis of the model. 

Let us also mention that, due to the absence of a spin operator in the charge density, the spin 
flip GP p( 01 ' 01 ) 1 would also be zero if we had neglected the energy difference between J = 1/2 
and J = 3/2 partners. This is why this GP is numerically so small with respect to its non spin 
flip partner. 

In the case of the p( n ' n )° polarizability we show separatly the contribution of the sum over 
excited states (referred to as para-magnetic) and the contribution of the seagull interaction (referred 
to as dia-magnetic). For the other GPs, the seagull does not contributes. 

7 Conclusion 

In this paper we have proposed a model independent analysis of low energy (e, e', 7) reactions. 
The output is concentrated in - a priori - 10 generalized polarizabilities which are functions of 
q, the virtual photon momentum. The only hypothesis actually needed for this analysis is the 
gauge invariance of the hadronic tensor. The decomposition of this tensor in two independently 
gauge invariant parts, one of which is exactly known (Hb) while the other one is regular (Hnb) 
is not unique. A regular gauge invariant term can be traded from Hb to H^b- This amounts to a 
redefinition of the polarizabilities and to avoid ambiguities or confusions we have been careful in 
defining Hb and Hnb- Our view is that while our decomposition is not unique it is well defined 
and does allow a meaningful comparison with any theoretical model. The reason is that, if the 
ten GPs are actually measured, then the complete VCS amplitude can be reconstructed up to and 
including order q'. 

To actually measure the GPs, we have proposed a method to analyze the threshold (e,e',7) 
reactions. The value of this analysis is that it handles exactly the interference between the VCS and 
the BH amplitudes, which is not simple because exchanged photon is not the same in both processes. 

To guide the analysis of the first experiments, we have made an estimate of the GPs in the 
non relativistic quark model. This is not a complete calculation since we have used simplifying 
approximations such as the neglect of recoil contributions. The latter increase with q and should 
be considered when q is larger than a few hundreds of MeV. 
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We also have omitted the contribution of the tt exchange in the t-channel. The reason is that 
in the NRQM, the quark mass is assumed to be large, thereby explicitly breaking chiral symmetry. 
So there is neither a compelling reason nor a clear way to introduce the pion coupling in this 
model. Moreover we recall that, according to duality arguments [17], the sum over s-channel 
resonances contains implicitly the t-channel exchanges. If an effective model like the NRQM bears 
some resemblance with reality, one may expect some double counting if one adds the t-channel 
exchanges in an uncontrolled manner. 

On the other hand we do know that the 7r° exchange is present in the experimental GPs and, 
due to the small mass of the pion, this may induce rapid variations as a function of q. In an ideal 
experiment this would cause no problem but in the real case it may help to add the 7r° exchange to 
Hb and subtract it from H^b- The 7r° exchange amplitude is by itself gauge invariant and regular 
in the soft photon limit [18]. Therefore this is the kind of term that one can freely trade from the 
Born to the non Born amplitude. This amounts to a well defined redefinition of the GPs and in 
this way the potentially dangerous variations due to the pion pole would be shifted to the known 
part of the amplitude, where they are innocuous. 

To conclude we do not wish to hide the fact that a fruitful analysis of threshold (e,e',7) 
reactions probably requires high quality data since the interesting information is obtained by a 
subtraction procedure, as explained in Section 5.1. This may be an interesting demonstration of the 
new electron accelerators which are now coming into operation. 
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Appendix A Conventions and units 

The units used in this paper are such that h = c = 1 and the magnitude of the electron charge 
is e = yj 47r olqeDi with olqed — 1/137. 

The metric tensor is {g^ v [i,v = 0,3) with signature (1,-1,-1,-1). Space components are 
labelled with Latin indices. 

Spherical harmonics, and Clebsch-Gordan coefficients follow the conventions of Edmonds [13]. 
The convention for the time reversed states is 

\k , jm > = (-) j+m \ -k,j-m> (133) 

for a particle of intrinsic spin j and rest frame projection m. 
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Appendix B Polarization states and spinors 



The unit vectors (e(i), i = 1, 2, 3) of the CM frame are defined in the text 
The polarization states of a photon with momentum k are 

k x e (1) 



e[k ,1) = e[k ,2) x k , e(k,2 



k x e (1) 



, e Mfe,3 



k, 



and the helicity states are as usual 

e(k,±i) = t 1 



e Ik, 1 1 ± i£ k,2 



£ lk,0] = £ f/fe,3 



° 1 *>(o) = ( k/ *° 



\/2 

The 4-vectors polarization states in the Lorentz gauge are then defined as 

£ "(±1) 
The helicity spinors are 



s/k Q + mxh[k 

2h^k - m Xh (k 
cos a/2 
sin a/2 



e-** sin a/2 
cos a/2 



where (a, <p) are the polar and azimuthal angles of the direction k. 
The rest frame spin-projection spinors are 

These spinors are the positive energy solutions of the Dirac equation (7.^ — m)u(k) 
Dirac matrices are those of ref.[19]. 
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Appendix C Vector basis 

The vector spherical harmonics are defined according to 

= E < U,WM > Yl(k) e (/x), (139) 
and the magnetic, electric and longitudinal vectors of the multipole expansion are, respectively 



-XLL ?L /""___ ~\)LL—1 , /_^_ 



- 1 










1 2L+ 1 




/ £ 


^ - 1 




V 2L+ 1 


/2L + 1 



LL+1 



(140) 



v -r x y " 
They have the following useful properties 

M L M (k) = L k Y&(k)/y/L{L + l), £% I = -ikxM L M (k), E L M {k) = kY^(k), (141) 

where L k is the angular momentum with respect to k . From eqs.(140, 141) it is clear that the 
electric and magnetic vectors are transverse and vanish if L = 0. 
The 4-dimensional basis is defined as 

V>( LM,i) = (^f) ) , = 

V"(lLM, S) = ( M ?^ ) , F"(2£M, J) = ) <142) 

and one has the closure relation 

9 pp V>1 (p LM > F "* (p lm > = 9^ 6 (k~ &') (143) 
and the orthogonality relation 

y dkV^ (pLM.k^jV^p L' M\Pj = g pp >8(L,L')8(M,M'). (144) 
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Appendix D 



In the non relativistic quark model defined in the text, the non zero GPs have the following 
expressions 

1 1 e -? 2 /6« 2 



p(0101)5 

p(0112)l 

(1111)5 
r vara 



dia 
p(01,l)5 

p(01,l)5 



(01,1)5 



rvS,Jx 

L d 



+ 



25+118 o? x= j^ A ^\ m ~ m X ' E{q)-E X {q) 



1 /3 e -« 2 /6« 2 



108 V 5 m.a 2 

9 X=jV*,A* 



V- 2 ("I) 

2_> a * 



7 2,S,Jx 
'ad 



r 7 2,S,Jx 



21 x 



m-mx E(q) - Ex(q) 



1 4 e-* 2 /6« 2 



7 5 



+ 



yS 



25+127 m 2 . \m-mA ' E(q)-E^(q) 



D (llll)5 . 7^e-5 2 / 6Q: 



54 m g ai 2 

(01,1)5 p (01,l)5 
^J 1 + -^5 

1 1 e-? 2 / 6 " 2 

~ 25 +1108 m Q a 2 



E 



2Z c 5,Jx 



X=JV*,A* 



) 



-g 2 /6a 2 



25 + 136^3 m q a 2 



E 



2 ("I) 



ix-i/2 ( 2j ljSjJx 

I ad 



X=N*,A* 



21 x 



m - mx E(q) - Ex{q) 



(145a) 
(145b) 

(145c) 

(145d) 
(145e) 
(145f) 

i(145g) 



where the P wave excited states (TV*, A*) taken into account are shown in the first row of table 
(2). The next rows give the coefficients ax taken from ref. [16]. The isospin of the state X is 
denoted Ix and the mass of the quarks and the oscillator parameter are, according to ref. [16], 
m q = 350 MeV , a = 320 MeV. The Z coefficients come from angular integration and spin sums. 
Their values are given in table (3). 



X 


7V*(1535) 


JV*(1650) 


JV*(1520) 


JV*(1700) 


A*(1620) 


A*(1700) 


ax 


0.85 


-0.53 


0.99 


0.11 


1.0 


1.0 



Table 2 Mixing coefficients taken from ref.[ll]. 

The seagull term contributes only to p( n ' n )°. The contribution of the A(1232) to this 
polarizability is referred to as para(magnetic), by analogy with the real Compton limit and the 
seagull contribution which has the opposite sign is referred to as dia(magnetic). 
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Table 3 Coefficients resulting from the angular integration and spin sums. 
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Table caption 



Table 1. Allowed quantum numbers of the non Born multipoles, to order q'. 
Table 2. Mixing coefficients taken from ref.[ll]. 

Table 3. Coefficients resulting from the angular integration and spin sums. 



Figure caption 



Figure 1. (a): FVCS amplitude, (b,c): BH amplitudes. 
Figure 2. Generalized polarizabilities in the NRQM. 



Final photon 


Initial photon 


Spin flip type 


P' = l 


p = 


L = 0,2 


5 = 1 


p = 2 


L = 2 


S = 1 


p=l 


L = 1 


5 = 0,1 


p' = 2 


p = 


L = 1 


S = 0,1 


p = 2 


L = 1 


5 = 0,1 


p=l 


L = 2 


5 = 1 



X 


JV*(1535) 


JV*(1650) 


JV*(1520) 


JV*(1700) 


A*(1620) 


A*(1700) 


ax 


0.85 


-0.53 


0.99 


0.11 


1.0 


1.0 



L 


S 


Jx 


d 




^L,S,Jx 
ad 


r^LjSjJx 


Z A 

^»1232 


1 





111 




a/273 


-2/^3 







1 





312 


2^/2/3 


2^/2/3 


2/^3 


-2/^3 




1 


1 


Ml 


2 


-2/3 


-2^2 


-2^2/3 





1 


1 


3/2 


-2 


2/3 


-V2 


-a/2/3 


1 


2 


1 


3/2 








^30 








0.00 



0.10 



-0.02 - 



0.00 



-0.04 



CL 

CD 



p (01,01)0 

(01,01)1 

100xP 



CL 
CD 



-0.10 



-0.06 



p(11.11)1 




p(11.H)0 


(para) 


□ p< 11 ' 11 >° 


(dia) 


p(11,11)0 


(total) 



-0.08 



0.0 



200.0 



400.0 600.0 

q (MeV) 



800.0 1000.0 



-0.20 



0.0 200.0 400.0 600.0 800.0 

q (MeV) 



1000.0 



0.0000 




-0.0030 ' ' 1 1 1 1 1 1 1 1 1 

0.0 200.0 400.0 600.0 800.0 1000.0 

q (MeV) 



0.0000 



-0.0005 

E 

CL 

CD 

-0.0010 



-0.0015 




0.0 200.0 400.0 600.0 800.0 1000.0 

q (MeV) 



